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We investigate the effects of pseudoscalar-photon mixing on electromagnetic radiation in the pres- 
ence of correlated extragalactic magnetic fields. We model the Universe as a collection of magnetic 
domains and study the propagation of radiation through them. This leads to correlations between 
Stokes parameters over large scales and consistently explains the observed large-scale alignment of 
quasar polarizations at different redshifts within the framework of the big bang model. 



I. INTRODUCTION 



Light pseudoscalar particles, such as the axion, arise naturally as pseudo-Goldstone bosons in spontaneously broken 
global symmetries [1-9 . An axion has an effective coupling to two photons, therefore in an external magnetic field it 
can oscillate into a photon and vice versa [10 21 . This mixing between axions and photons can lead to observable 
changes in the intensity and polarization of photons, having interesting astrophysical and cosmological implications 
[22TI37] . Since the mixing depends on the frequency of radiation, it also affects the electromagnetic spectrum [38l - 
21] • Various experimental searches of pseudoscalar particles have led to significant limits on their masses and on the 
coupling parameters of the models CHI 1231 H5H6T] . 

In the current paper we study pseudoscalar-photon mixing in the presence of correlated background magnetic fields 
in order to understand the results reported by Hutsemekers et al. [62H64] pertaining to the coherent alignment of 
quasar polarizations over Gpc scales. It has been shown earlier that propagation of radiation through the extragalactic 
medium, in the presence of a background magnetic field, can affect all Stokes parameters due to pseudoscalar-photon 
mixing. These changes have been investigated for radio [HI E3 ES] , optical [JTJ I52T[5H loTHST)] , and cosmic microwave 
background (CMB) [33l 021 EHl [70l E] photons. We explicitly show that in our model it is possible to obtain nonzero 
correlations between the optical polarization of quasars separated by large distances in the sky. These nonzero 
correlations lead to an effect similar to the observed large-scale alignment of the polarizations. 

We consider the extragalactic medium to be a large collection of magnetic domains, the background magnetic field 
in each domain being constant. This model was recently used in [71] where the extragalactic medium was considered 
as a large collection of uncorrelated magnetic domains in order to study the effects of pseudoscalar-photon mixing on 
CMB polarization. We extend the model of [7T] to include correlations between magnetic fields in different domains, 
as motivated in |72H74j . and study the propagation of optical radiation from quasars through these domains. 

The paper is organized as follows. In Sec. [II] we outline the main concepts and equations of pseudoscalar-photon 
mixing and describe the domain propagation model. In Sec. |III|we obtain magnetic field correlations between different 



domains. In Sec. IV we bring together results derived in the previous two sections and study alignments among the 
angles of polarization using a full numerical propagation model. We conclude with a brief discussion of our results and 
offer perspectives in Sec. [V] We also discuss an approximate analytical treatment to explain correlations in quasar 
polarizations in the Appendix. 

II. PSEUDOSCALAR-PHOTON MIXING 

In this section, we consider the coupling of a light pseudoscalar to an electromagnetic field and discuss the prop- 
agation of electromagnetic waves in the presence of a background magnetic field. We first present a short review of 
the field of pseudoscalar-photon mixing in |II A| and subsequently describe our model of propagation in |IIB| 

A. Basic concepts and equations 

The interaction Lagrangian for the coupling of pseudoscalars with an electromagnetic field can be written as, 

£ int = ^F^F^, (1) 
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where g<f, is the pseudoscalar-photon coupling constant, <\> the pseudoscalar field, F^ v the electromagnetic field tensor, 
and F^v = \^^, V pijF f " T its dual. A single pseudoscalar is therefore effectively coupled to two photons. As shown in 
[TBI IT5] , the longitudinal component of the background magnetic field plays a negligible role in pseudoscalar-photon 
mixing. Since only photons polarized parallel to the transverse component of the background magnetic field (St) 
decay this effect can lead to a change in polarization of the electromagnetic wave. Alternately photons can mix with 
off-shell axions or axions may decay into photons, in either case leading to a change in polarization. Pseudoscalars can 
thus mix with photons from distant galaxies or the background CMB radiation, and lead to a rotation of polarization. 

In order to describe the propagation of electromagnetic waves, we choose a coordinate system such that the z axis 
lies along the direction of propagation, and the x axis is parallel to the transverse component of the background 
magnetic field Bt- We define the gauge invariant quantity A — E/lu (E being the usual electric field and ui the 
frequency of radiation) and resolve it in components parallel (^4y) and perpendicular (A±) to the direction of Bt- 
This enables us to write the field equations solely in terms of Ah and <f>, since A±_ does not mix with the field <j), 



+ 9 l) lu\ - M = °- (2) 



Hz) 

The "mass matrix" or "mixing matrix" , M, is given by, 

m=( u i -^y 0) 

where uip is the plasma frequency, the pseudoscalar mass, and Bt = \Bt\- 

The solution of the field equations follows from diagonalizing the above matrix equation using the mixing angle 9 
defined by, 

tan20 = / 5<A S T , (4) 

where I denotes the oscillation length, 

l=- 2 g-. (5) 

up -m% 

In this paper we will be working in the limit of very small pseudoscalar mass, <C up, since for masses much heavier 
than this, the mixing with photons produces a negligible effect for intergalactic propagation, for the range of allowed 
parameters. 



B. Equations of polarization propagation 

Polarization of radiation can be quantified in terms of Stokes parameters, which can be written as linear combi- 
nations of correlations between different components of the field A. A convenient way to formulate the correlation 
functions between initial components of A and <j) is to arrange them as elements of a physical density matrix p(0), 
given by, 

/ (4,(0)^(0)) <A||(0)A1(0)> (A||(0)^(0)) \ 
p(0)= <4i(0)Afi(0)) (Ax(O)^l(O)) (Al(O)0*(O)> . (6) 
V (^(0)^,(0)) (0(0)^1(0)) (0(0)0* (0)) / 

The correlation functions propagated through a distance z can then be expressed as, 

P (z) = P(z)p(0)P(z)-\ (7) 

where the unitary matrix P(z) describes the solution to the field equations, for a given mode uj |18|, I71j . In our 
coordinate system, P(z) is given by, 

(1 — 7 sin 2 9 7 cos 9 sin 9 \ 

e - l [c+A A -(^-^) 1 /2] z , (8) 

7 cos 9 sin 9 1 — 7 cos 2 9 J 
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where 7 = (1 — e' lAz ), A — A^ — Aa, and Aa, A^ are defined in terms of the frequency, u, and the eigenvalues, p,±, 
of the matrix M, 

A A = ^Juj 2 - (9) 
A 4 , = Ju?^K-w. (10) 

For typical values of the electron density (n e m 10~ 8 cm ), and for optical radiation of quasars, with frequency 
v w 10 6 GHz (which corresponds to I w 4 Mpc, here v = uj/2it), we have ui 3> wp, m^, and g^Br- Also, for typical 
values of the coupling, = 6 x 10~ n GeV -1 , and magnetic field, Bt = 1 nG, we can approximate A as, 

A = A -A j4 wyv / 1 + tan 2 29 = - sec 26. (11) 

We now study the propagation of electromagnetic radiation through the intergalactic medium in the presence of 
pseudoscalar-photon mixing. Similar studies have been carried out earlier by several authors, e.g. [28 30, 33, 36, 71j. 
It is reasonable to model the medium as a large number of magnetic domains, with a uniform direction and strength 
of magnetic field in each domain. Additionally, the medium is assumed to have a uniform value of uip. In the ensuing 
analysis we will use this model to explain correlations in the optical polarization of quasars as reported by Hutsemekers 
et al. |62H64| . A crucial component of our analysis will be the inclusion of correlations between magnetic fields in 
successive domains which has hitherto not been considered in any such study to the best of our knowledge. Details 
of these correlations will be presented in the next section. 

The transverse magnetic field TBt in the i th domain is taken to be oriented at an angle (3i with respect to the x 
axis (the "parallel axis") of the external coordinate system. Starting with the density matrix p(0) = diag(l,0,0), 
corresponding to an initially unpolarized electromagnetic wave, the density matrix is propagated through each domain 
using |7]) . After propagation through each domain, the electromagnetic wave vector is rotated back in order to account 
for the change in direction of the transverse magnetic field from one domain to another. The propagation through n 
domains of size z each (the n th one being closest to the Earth) gives us an expression for p n (Z = nz) |71j . 

Pn {Z = nz) = R- 1 {p n )P{z)R{p n )R- 1 {p n ^ 1 )P{z)R(p n ^ 1 )R- 1 {p n ^ 2 ) ... R- 1 ((5 1 )P(z)R(l3 1 ) 

x p(0)R- 1 ((3 1 )P- 1 (z)R((3 1 )R- 1 ((3 2 )P- 1 (z)R(/3 2 ) ... R- 1 ^ n )p- 1 (z)R((3 n ), (12) 

where R(f3 m ) is the rotation matrix that acts only on the two-dimensional space transverse to the propagation 
direction; i.e., it represents a rotation by the angle /3 m about the z axis. Explicit expressions for the propagation are 
given in the appendixes of [71] . 

We also give here expressions for the reduced Stokes parameters in terms of different components of the density 



matrix, for use later in the paper, 

I(z) = p 11 (z)+p 22 (z) = {A ll (z)A^(z)} + (A ± (z)Al(z)}, (13a) 

Q(z) = p n (z)-p 22 (z) = (A ]] (z)AJ l (z))-(A ± (z)Al(z)), (13b) 

U(z) = p\ 2 {z) + p 2 i(z) = {A\\(z)A* L {z)) + (A±(z)A*\(z)), (13c) 

V{z) = i( Pl2 (z)-p 21 (z)) = i((A ll (z)A* x (z))-(A ± (z)A^(z))). (13d) 

Also, the linear polarization angle if) and the degree of polarization p are given in terms of Stokes parameters by, 

tan2^ = U/Q, (14) 

p = VQ 2 + U 2 + V 2 /I. (15) 



III. MAGNETIC FIELD CORRELATIONS 

We now calculate the correlations between components of the magnetic field in different domains. The analysis 
derives from the magnetic field spectrum M(k) discussed in (72HH], which can be defined using magnetic correlations 
of the form, 



(bi{k)b*(q)) = dk, q Pij{k)M{k) 



(16) 
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where crfj(k) = Py(fe)M(fc) and k = \k\. Here fo(fe) is the Fourier transform of the present day magnetic field B(r) 
and Pjj(fe) = — ^jpr^ is the projection operator included to be consistent with a divergenceless magnetic field. 
Although this distribution was originally proposed for primordial magnetic fields, it is justified to use it here since on 
galactic and larger scales the magnetic field simply redshifts away, as discussed in [75H76"] . The magnetic field in real 
space is defined as, 



V 

where V is the volume in real space. In accordance with [72], we consider a power-law dependence for M(k), 

M(k) = Ak nB ; n B > -3, 
where n B is the power spectral index, and use a sharp fc-space filter (window function) of the form, 



W(6) 



9 < 1 
9 > 1 



(17) 



(18) 



(19) 



Using the above relations and taking the continuum limit by replacing ^2 k by / d 3 k, we find that the spatial 
correlation between magnetic field components at two different points in space separated by a distance r' is given by, 



(B i (r + r , )B i (r)) 



1 

V 



d 3 k 
(2tt)3 



e ik - r 'aUk)W*(kr G ), 



(20) 



where tq is the "galactic" scale, taken to be 1 Mpc here. The constant A in ( 18 1 is evaluated by assuming a smooth 

„-i i ». „-i 



variation of the field B(r) over the scale kc 



1 Mpc" 



Further, we consider a completely isotropic distribution 
of magnetic fields for simplicity, for which Pij(k) = (2/3)<5y. In view of the above assumptions, the spatial correlation 
calculated over a sphere of radius tg = 1 Mpc gives the constant A as, 



Vn 2 B 2 {3 



J.3+TII 



(21) 



where Bq is the fiducial constant magnetic field whose value is assumed to be 1 nG [741 177) . It may be appropriate 
to impose a large distance cutoff, r max , on the correlation for distances comparable to the size of the system (the 
Universe). In such a case the lower limit of the integral in (20 1 will be £; m ; n = r~^; here we assume that this cutoff 
is small enoug h an d can be approximated as zero. 

We now use ( 20 ) to obtain the correlations between magnetic field values in different domains. Consider the m th and 
n th domains along a given line of sight, separated by the distance r mn . The transverse magnetic fields (of strength Bt) 
in these domains are aligned at angles /3 m and j3 n with respect to the fixed external coordinate system, respectively. 
Hence the correlation between their components in the fixed coordinate system can be computed as, 



B 2 T - 



cos/3 

'ill 

cos/3 



dk 



sin kr TI 



(22) 



using Pij(k) — 5ij (since for transverse magnetic fields fcj 



0). Here Ag(fc) 



n B 



(4) 



3+71 E 



J__k 
V 2 



^M(fc) = f (3 + 



is the power per logarithmic interval in fc-space residing in the magnetic field. Therefore, 



B% (cos fim COS (3 n 



-B 2 
2 \r. 



TG 



3+rif 



(3 + ns) 



A+n B sma , fix 



(23) 



where the term in square brackets is of order unity. For example, for ng = —2.37, which was reported as the best fit 
value in the statistical and numerical analysis of the matter and CMB power spectrum in [77] , an explicit numerical 
evaluation yields a value for this term between 0.96 and 1.33 for r mn /rG — 1 and large r mn /rG, respectively. Further, 
for ns — > —3.0 this term rapidly approaches unity. 

For the completely isotropic case considered in this analysis we have -Py(fe) oc Sij, leading to a cancellation of any 
cross correlations between transverse components of the magnetic field, i.e., (cos/3 m sin/3 n ) = 0. Also, the invariance of 
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correlations under rotation of the external coordinate system immediately implies an equality between (cos /3 m cos /3 n ) 
and (sin/3 m sin/3„). The final set of correlations can therefore be summarized as: 

1 B 2 ( \ 3+ ™ B 

(cos m cos /U = dr — , (24a) 

(cos/? m sin/?„) =0, (24b) 

1 d2 / . \ 3+ri B 

(sin/? m sin/U = ^dr — ■ ( 24c ) 



IV. CORRELATIONS IN THE OPTICAL POLARIZATION OF QUASARS 



In this section we use the domain propagation model described in Sec. [TTJ with magnetic fields correlated according 
to the correlation functions discussed in Sec. |III| to understand the large-scale alignment of quasar polarizations 
observed by Hutsemekers et al. [52Tf6"i] . We perform a detailed numerical analysis of the correlations by direct 
integration of the propagation equations, using a magnetic field distribution that follows from our discussion in the 
previous section. We also refer interested readers to the Appendix, where we have included an analytical approach to 
estimate such correlations in quasar polarizations in the regime z/l -C 1 and 9 <C 1. 



We first extract a one-dimensional (ID) magnetic field distribution along a given line of sight from (20), by inte- 
grating over the transverse components of k. This gives 



(Biiz + ^Bjiz)) = 



ka 



dk z 
2tt 



where ajAk z ) = SijM(k z ). The ID magnetic field spectrum M(k z ) is now given by, 



M{k z ) 



A 



,2+n E 



3ir \ 2 + riB 

A distribution function compatible with the above correlations can be written as, 



•&?(*« 



2M(k z ) 



(25) 



(26) 



(27) 



where Ni(k z ), Nj{k z ) are the normalizations. This represents an uncorrelated Gaussian distribution for two compo- 
nents of the magnetic field in Fourier space, corresponding to the wave vector k z . 

Using the above distribution, we generate the transverse Fourier components b x (k z ) and b y (k z ) of the magnetic 
field for the discretized wave vector k z . This can be done for each domain in Fourier space independently as the 
distribution is uncorrelated for different k z . The magnetic field components B x {z) and B y {z) in real space can then 
be obtained by implementing the inverse Fourier transform in (17). The transverse component of the magnetic field is 

characterized by the magnitude Bt = B 2 (z) + B 2 (z) and the angle j3 — tan -1 (B y (z) / 'B x (z)) that it makes with the 

x axis of the external coordinate system. We repeat this exercise for every domain to obtain the complete magnetic 
field landscape along a given line of sight. 

We point out that the magnetic field changes relatively slowly over small distances of order 10 Mpc. This can be 
seen in Fig. [l] where we plot a random realization of one of the components of the correlated magnetic field. Hence 
our line of sight calculation applies approximately even if quasars lie at slightly different angles. For larger angular 
separations we require a full 3D simulation which is very time-consuming and we postpone this for future research. 

Next, using the above model, we study propagation along a given line of sight through 2 12 domains of size z = 1 
Mpc each, with 400 equally spaced quasars - the nearest one being 10 Mpc away and the farthest one being 4 Gpc 
away from the Earth. The initial electromagnetic radiation of each quasar is assumed to be unpolarized and the 
parameters and v are chosen to have typical values 6 x 10~ n GcV -1 and 10 6 GHz, respectively. The simplest 
statistic to test whether our model predicts an alignment of quasar polarization vectors is to calculate the dispersion of 
the polarization vector of a quasar with respect to its nearest neighbors [67] . For this exercise, we group quasars into 
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FIG. 1: A random realization of one of the components of the correlated magnetic field as a function of distance. 



l) , 



sets of n v quasars each, and calculate the mean value of the vector [cos(2-0j), sin(2^)], ipi being the polarization 



angle [calculated using (14)], for each set. A measure of the dispersion d k of vectors in each set is given by, 

d fe = — Vcosp^-Vfc)], (28) 
i—i 

where is the mean polarization angle in the k th set. The statistic defined as, 

1 



S£ = -}^d k , (29) 



gives us a measure of the dispersion in the data. A large value of S^, indicates a strong alignment between polarization 
vectors. 

In Fig. [2]we show results of the S£ statistic with n v = 10, obtained using polarization angles from our propagation 
model and compare it with data generated using a random distribution of the polarization angles. Also, we have 
checked that the distribution of polarization angles obtained using random magnetic fields, when quasars along a 
given line of sight are at slightly different angles so that each experiences completely random magnetic fields on 



propagation, mimics a random distribution of polarization angles. Since the mean value of the statistic in (29) 
is clearly larger for data obtained using the propagation model, this provides conclusive evidence of a preferred 
alignment among quasar polarization vectors. Further, we expect the mean value of this statistic for a random sample 
to be proportional to l/^/n^ [57]. We indeed find that results for the random sample are well described by the fit 

S£ « 0.90/ \fn v (Fig. |3j). Further, for the propagation model, S£ ~ 1.26/n° 47 . This approximately behavior 
in our model suggests that the angles of polarization are preferentially aligned in each neighborhood, yet span the 
whole range of possible angles as in the case of a random distribution of angles. 

A related effect of correlated magnetic fields in our propagation model is that they lead to a larger axion conversion 
probability, defined as 



(<f>{z)<p{z)), (30) 



compared to that obtained using a model with random magnetic fields. Therefore it is more likely for photons to 
decay into axions, over a given distance, in the presence of a correlated magnetic field distribution. This can be seen in 
Fig. [4] where we show the development of the conversion probability with distance for an initially unpolarized quasar, 
100 Mpc away from us, as optical radiation travels through correlated and random magnetic field distributions. 

We finally discuss how our results change with a different choice of kc or equivalently tq. Here we assume 
that the magnetic field is uniform in each domain, and the size of each domain is equal to tq. In the results 
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FIG. 2: Histograms of the statistic So with the number of nearest neighbors n v = 10, for 10 4 random samples (left-red) and 
10 4 realizations of the magnetic field distribution in the propagation model proposed in the paper (right-blue), using parameter 

values n e — 10~ 8 cm -3 , nj = —2.37, v = 10 6 GHz, and g$ = 6 x 10" 11 GeV" 1 . The larger value of the mean 5*^ obtained 
for the propagation model, compared to the random sample, indicates preferential alignment of quasar polarization vectors, as 
reported by Hutsemekers et al. [62-64 . 
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FIG. 3: The mean (left panel) and standard deviation (right panel) of the statistic Sp as a function of nearest neighbors n v 
for 10 4 random samples (red circles) and 10 4 realizations of the magnetic field distribution in our model of propagation (blue 
squares), using identical parameter values as Fig. [2] The variation of the mean S£ with n v has been plotted on a log- log scale. 
The lines are fits generated using a power-law dependence of the form ax b [where b « —0.50, a w 0.90 for the random sample 
(red); b w —0.47, a « 1.26 for the propagation model using correlated magnetic fields (blue)]. 



presented so far we have chosen tg — 1 Mpc. A larger value of r G will induce correlations over larger distances, 
hence we would naively expect a stronger signal for larger r G or equivalently smaller k G . However, the phenomenon 
is nonlinear and complicated. For example, the photon to pseudoscalar conversion probability does not increase 
monotonically but oscillates. This is because of the reverse process of pseudoscalar to photon conversion, once the 
population of pseudoscalars becomes sufficiently large. Furthermore the orientation of linear polarization does not 
evolve proportionately to the conversion probability. Hence it is not possible to anticipate how the results depend 
on ka- In Fig. [4] we also show the conversion probability for k G = 2 Mpc -1 and k G = 0.5 Mpc" 1 . We see that 
the conversion probability slightly decreases in both cases. On calculating the statistic for these cases we find 
that for n v = 10 nearest neighbors, S£ w 0.43 for k G = 1 Mpc" 1 , 0.44 for k G = 2 Mpc" 1 , and 0.48 for k G = 0.5 
Mpc -1 . Therefore for larger r G there are higher correlations although the conversion probability is smaller due to the 
conversion of pseudoscalars into photons. For smaller r G , on the other hand, the correlation is only feebly affected, 
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FIG. 4: Axion conversion probability as a function of distance traveled, for optical radiation from a single initially unpolarized 
quasar 100 Mpc away from us, as it encounters correlated (solid thick black curve for kc = 1 Mpc -1 , small-dashed black curve 
for kc = 0.5 Mpc -1 , large-dashed black curve for ka = 2 Mpc -1 ) or random (solid thin black curve) magnetic field distributions 
on its journey. Each curve is generated with parameter values n e — 10 -8 cm -3 , ub — — 2.37, v = 10 6 GHz, g<t> = 6 x 10 -11 
GeV -1 , and is averaged over 10 4 realizations of the respective magnetic field distribution. 

and the conversion probability decreases for most values of the distance, as expected. 

V. DISCUSSION 

A very striking alignment in the optical polarization of quasars has been observed over cosmologically large distances 
(~ 1 Gpc) by Hutsemekers et al. [52146~4"] at both low (z ~ 0.5) and high redshifts (z ~ 1.5). It is unlikely that quasar 
polarizations are intrinsically aligned with one another over such large distances since this appears to be in conflict 
with the basic assumptions of isotropy and homogeneity of the Universe. It has been proposed in the past that these 
observations may be explained in terms of a propagation effect related to axion-photon mixing. In the current paper 
we have explored this possibility in detail, and studied the propagation of electromagnetic radiation through a large 
number of magnetic domains. 

We have shown that correlations between components of the background magnetic field in various domains induce 
correlations between the linear polarization of different quasars. Our analysis explicitly shows that quasar polarizations 
are expected to be aligned in such a scenario due to the presence of significant correlations between them, as observed 
by Hutsemekers et al. [62 -64 . A calculation of the dispersion of polarization with respect to nearest neighbors in a 
model of propagation along a given line of sight shows significant alignments compared to a random sample. A direct 
comparison with observations and further study of the effect clearly requires numerical simulations of the complete 
3D propagation model that we have proposed. We leave such an analysis to future work. 

We further note that a calculation of correlations in the circular polarization, Stokes V parameter, is also expected 
to yield a significant result. A detailed calculation of this parameter is postponed to future research. This may prove 
to be an important signature of pseudoscalar-photon mixing in the observed alignment of quasar polarizations, as 
noted in [75] , 
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Appendix: Analytical treatment of the correlations 

In this appendix we discuss an analytical approach to estimate the correlations in quasar polarizations for small 
domain size (zjl <C 1) and small mixing angle (0 <C 1). In this limit, we expand the 22 element of QSJ) to leading 
order, 

e -i[ u+ A A -(u, 2 -u,%y/*]z _ j _ f * ^ _ sec2 6»), (Al) 

where we have omitted the overall phase factor. We can further write the propagation matrix as the sum of a quasiunit 
matrix and a matrix proportional to sin0, 

P(z)=l{z)+$m6V{z,6). (A2) 
Suppressing the irrelevant overall phase in Q, the matrices I(z) and V(z,9) are given by, 

Z(z)=diag(l,l,e lAz ), (A3) 

and, 

(—7 sin 7 cos \ 

if sin0sec20 . (A4) 
7 cos0 7 sin0 / 

We can now express p n {Z = nz) as a power series in sin0. For an initially unpolarized electromagnetic wave 
corresponding to a density matrix p(0) = diag(l, 1, 0), the calculation outlined above gives the following leading order 
terms in different powers of sin 9 V(z,8), 

pM(Z) = p(0), (A5) 

n 

p%\Z) = 5>M(Z,i), (A6) 



n k—1 



P [ n\Z) = Y^zZpn^j^ ( A7 ) 

k=2 j=l 



where, 



/ -2 sin 2 9 cos 2 ^(1 - cosAz) - sin 2 9 sin 2^(1 - cosAz) 7 * e - iA *( n -J') cos/3,- sin(20)/2 \ 
pW(Z,j)= - sin 2 0sin2/3,(l- c °s Az) -2 sin 2 sin 2 (3 3 (1 - cos Az) 7 * e - iA *(*-i) sin/3,- sin(20)/2 , (A8) 
\ 7 e iA ^"-J) cos ft- sin(20)/2 7 e iAz ("-J') sin/3 3 sin(20)/2 / 

and, 

Pn.ii i' fc ) = _4 sin2 ^ cos2 ^C 1 _ cos Az ) cos[Az(fc - j)} cos /3 fc cos ^ , (A9a) 
p l „\ 2 (Z, j, k) = -2 sin 2 cos 2 0(1 - cos Az) cos[Az(fc - j)} ( sm(f3 k + (3 3 ) - i sin(/3 fc - fa) tan[Az(fc - j)]) , (A9b) 

Pn] 21 (Z, j, k) = -2 sin 2 cos 2 0(1 - cos Az) cos[Az(fc - j)] ( sin(/3 fc + /3j) + i sin(/3 fc - ^) tan[Az(fc - j)]J , (A9c) 

p [ ^ 22 (Z,j,k) = -4sin 2 0cos 2 0(l - cos Az) cos[Az(/c - j)] sinft sin/3j. (A9d) 

Here we only report elements of the 2x2 submatrix relevant to polarization for p n (Z,j,k). Further, we ignore 
higher-order terms in writing (A9). This approximation is valid as long as sin 2 and (z//)sin 2 are much smaller 
than unity. 

For two quasars (along a given line of sight) at distances Z± = nz and Z 2 = mz, where n and m are the number of 
domains of size z that lie between us and the two quasars, we wish to estimate the correlation (U{Z\ = nz)U(Z 2 = 
mz)). Using (13) we find that, 

(U{Z 1 )U{Z 2 )) = {PnM Z l)Pm,12( Z 2)) + (Pn,12(^l)pm,2l(^ 2 )) 

+ (Pn,2l(Zl)p m , 12(-^2)) + (Pn,2l(Zl)p m ,2l(Z 2 )). (A10) 
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On removing the constant zeroth order part of the propagated density matrix we can write, 

Pn (Z) a p m(Z) + pW(Z). 



(All) 



Using ( A8) and (A9) in this expression to calculate (AlOl, we end up with four-point correlation functions in the sines 



and cosines of rotation angles. We assume that the magnetic field fluctuations at any point are Gaussian, so that we 
can express these four-point correlations in terms of two-point correlations using the following decomposition, 



(S( ri )S(r 2 )S(r 3 )S(r 4 )) = (S( ri )S(r 2 ))(S(r 3 )S( n )) + {S(n)S(r 3 )) <S(r 2 )S(r 4 )) 4 

where S(r) is some function (such as a Gaussian) that falls off sufficiently fast with 
and write the correlation in (AlOl as, 



(S( ri )S(r 4 )){S(r 2 )S(r 3 )),(A12) 
r. Now we can further use (|24l. 



(fMZj ,r ( Z 2 )} 4{ ) sh. l 0(] -cosA:i- x 



n rn / \ 2(3+n B ) 

r G 



EE 



j n rn q—1 

2-Heee 



3+n [ 



q=2p=l \ r l* r W< 
n j — 1 rn q—1 



COS 



'EEEE 

j=2 i=l q=2 p=l 




i=lp=l 

rn n J — 1 

cos[A z (<? - p)] + E E E 



p=l i= 2i=l v r JfV 



3+n f 



cos[Az(j 



3+ni 



+ 



qi' pj 



3+rif 



( cos[Az(j + g — i — p)] + cos[Az(j — q — i + p)]j 

(A13) 



The above equation suggests that the linear polarization of quasars at different redshifts along a given line of sight 
is correlated, as observed by Hutsemekers et al. [62H64] . Here in using the correlations (24) we have set the term in 
square brackets in (23) to unity. Further we choose a constant value of 1 nG for Bt- These approximations affect the 



analytical result to within a factor of order unity. 

The analytical calculation presented here makes several approximations and hence only provides a rough estimate of 
the correlations. Nonetheless, it is very important since it establishes the presence of correlations among polarizations 
of electromagnetic radiation from quasars separated by large distances. These correlations arise due to the presence 
of correlations in the intergalactic magnetic field. The analytical result also provides useful insight into how the 
correlations among polarizations depend on the intergalactic magnetic field. In order to study correlations beyond the 
regime of validity of the analytical method discussed here, a complete numerical treatment of propagation, as done 
in Sec. |IV| of this paper, is necessary. 

It is also interesting to use the numerical propagation model in the regime where the analytical calculation holds, 
and compare the result with that obtained using (A13). For this, we study propagation along a given line of sight 
through a total of 2 10 domains of size z = 1 Mpc each, with 100 equally spaced quasars - the nearest one being 
10 Mpc away and the farthest one being 1 Gpc away from the Earth. The initial electromagnetic radiation of each 
quasar is assumed to be unpolarized. We choose the following parameter values: electron density n e — 10 -8 cm -3 , 
spectral index tib = —2.37, frequency of radiation v — 10 6 GHz, coupling = 6 x 10~ 12 GcV -1 , and Bo = 1 nG. 
This choice of parameters ensures that z/l < 1 and 9 « 1. Further, in the analytical result we set Bt = 1 nG and 



calculate as ri p = z\i — p\, while 



tq. In the numerical propagation, 



the value of Bt and the angle /3 in 
IV For the purpose of comparison 



each domain are obtained from the magnetic field distribution discussed in Sec. 
with the analytical result, however, we set the magnitude of the magnetic field at every point to be equal to 1 nG in 
the numerical calculation as well. Direct calculation shows that allowing the magnitude to vary makes a significant 
difference in the final answer. Hence it does not provide a fair comparison with the analytical result where we have 
fixed the magnitude of the field in order to make the calculation tractable. Fig. [5] shows a plot of analytical and 
numerical values of the correlation between the Stokes U parameters of two quasars separated by a relative distance 

C- 

We find that the analytical and numerical results agree to within about 20%. The discrepancy between the two 
results may be attributed to several other approximations made in the analytical treatment, where (i) we have set 
the term in square brackets in ( 23 ) to unity and (ii) we have ignored higher-order corrections. As a rough estimate, 
including the term in square brackets in ( 23 ) increases the analytical values of correlations by 60% — 70%. In addition, 
we also observe that the numerical propagation model has a small dependence on the "system size" (i.e., the total 
number of domains chosen), unlike the analytical result in (A13). This can be ascribed to dependence on the number 
of domains, of the discrete values of the wave vector k z used in the inverse Fourier transform. As the magnetic field 
values depend on k z , this introduces a small system size dependence in the results. Hence we should expect the two 
results to agree only within a factor of about 2. We also point out that we expect best agreement for small distances 
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FIG. 5: Comparison of analytical and numerical results for parameter values n e — 10 8 cm 3 , nj = — 2.37, v = 10 6 GHz, and 
= 6 x 10" 12 GeV" 1 . The plot shows the correlation {u(Zi = 100 Mpc)u(Z 2 = (100 + <) Mpc)) for different values of the 
relative distance £ between two quasars, where u is the normalized value of the Stokes parameter U with respect to the initial 
intensity, 1(0) = 2. Red circles show values obtained using the analytical result (A13l and blue squares are values obtained 
after the complete numerical propagation, averaged over 10 4 realizations of the magnetic field distribution. 



of propagation. This is because for larger distances the higher-order terms, dropped in the analytical calculation, 
cannot be neglected. In this case we will have to sum over a large number of domains and the higher-order terms will 
become increasingly important. The approximation used in (Al) will therefore break down. 
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